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SUMMARY 

By the method of superposition of conical flows, the load distribur- 
tion is calculated for regions of a long, rectilinear, svept-back wing 
behind the points at which the Mach lines from the trailing— edge apex 
intersect the leading edge. It is found that a good approximation to 
the load can be obtained by the application of a fairly simple correc- 
tion factor to the two-dimensional subsonic distribution. The si milar ity 
to two-di men sional flow is used to derive expressions for the loss of 
lift behind the Mach lines from the tips of the wing. 


INTRODUCTION 

In a previous paper (reference l), formulas were derived by the 
superposition of linearized conical flows (reference 2) for the loading 
of a thin, flat, swept— back wing of rectilinear plan form traveling at 
supersonic speeds, under the limitation that the Mach lines from the 
trailing-edge apex did not cut the leading edge. Ahead of these Mach 
lines the flow is known to follow the conical pattern identified with the 
triangular wing having the same apex angle. By superposing an appropriate 
distribution of other conical flows behind the trailing edge to cancel the 
lift in the wake, the triangular-wing loading can be corrected within the 
region behind the trailing-edge Mach lines to conform to the Kutta condi- 
tion for a subsonic trailing edge. This method is described more fully 
in reference 1. 

When the region of influence of the trailing edge includes a part of 
the leading edge and the stream ahead of it, the mutual interaction of the 
two regions results in a flow which is susceptible to treatment by the 
conical-flows method only at the expenditure of considerable effort. The 
flow in the tip regions is particularly difficult to calculate. A method 
of cancellation of external lift by supersonic doublet distributions, pro- 
posed in reference 3, offers some promise of providing a solution in the 
tip region if the pressure at all points inboard of the tip Mach cone is 
known from other sources. The method appears, however, to be in error as 
applied in the neighborhood of the leading edge or of a raked— out tip. 
Formulas of considerable simplicity have been derived by Heaslet, Lomax, 
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and Spreiter (reference 4'), but only for a particular family of slightly 
tapered wings with a filleted trailing edge and plan forms very slender 
relative to the Mach cone. 

Since no satisfactory method appeared to be available for investi- 
gating the general problem of the rectilinear swept-back wing with 
interacting leading and trailing edges — that is, wings of high sweep 
or aspect ratio, or wings of moderate sweep and aspect ratio at low 
supersonic Mach numbers — the work of reference 1 was extended, with 
some simplifying assumptions for the tip regions, to cover such cases. 

The principal object of the investigation was to obtain, if possible, an 
indication of the manner in which the flow in the outboard regions of a 
long narrow wing approaches the subsonic two-dimensional type of flow 
indicated by simple sweep theory. 

The calculations showed that the section loading assumes the general 
shape of the subsonic two-dimensional loading at the very station at 
which the trailing-edge Mach line intersects the leading edge. Further 
investigation showed that everywhere behind the Mach lines from the point 
of intersection and farther outboard the load could be closely approxi- 
mated by a fairly simple adaptation of the two-dimensional formulas. 

The following material outlines the method of obtaining the load 
distribution in the region of leading— and trailing-edge interaction by 
the superposition of conical flows and presents the formulas for approxi- 
mating the results by a single step based on correction of the two- 
dimensional formulas. A corresponding approximate correction to the load 
near the tips is also derived. 


CALCULATION OF LOADING BY SUPERPOSITION OF CONICAL FLOWS 

Procedure 

Figure 1 shows the central portion of a long,- narrow, swept-back 
wing divided into various regions by the trailing-edge Mach lines, the 
Mach line arising at the intersection of the trailing-edge Mach line 
with the leading edge, and so on. In the superposition process to be 
described, these Mach lines mark the regions of influence of the flow 
fields used in successive cancellations of lifting pressure outside the 
wing boundaries. 


. In region I, the streamwise component u of the upper— surface per- 
turbation velocity — to which, in linear theory, the loading is propo]>- 
tional — is given, as previously mentioned, by the known solution for the 
triangular wing swept behind the Mach cone 


UA(a) 


_mu c _ 

J m 2 -* 2 


( 1 ) 
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where u 0 , the constant velocity along the center line a=0, is 


u 


o 


mVa 

3 E( */l-m2) 


( 2 ) 


(See Appendix A for explanation of the symbols.) 

The flow over the remainder of the wing, however, is affected by 
the reduction of the lift to zero behind the trailing edge. The primary 
induction effects on the wing are given in reference 1 in the form of 
two corrections to the basic perturbation velocity: a symmetrical 

correction 1 


(Au)o , EMd 

0 0 K( 


with 


cp = sin 


-1 A-v 


V l-m^ 

and the integral with respect to a of 

(Au) = - ^ cos -i (_ 1 ~ a )( t a~ m t) ~ (mt-a)(l-t a ) 

(l-m fc )(t a -a) 


(3) 


(b) 


where u a = (du^/da)da. 


These induced velocities appear on the wing surface (in regions II, 
IIIj • • • ) as corrections to the basic velocity ua* Beyond the actual 
surface of the wing they represent a discontinuity in pressure which 
cannot be supported in the free stream and must therefore be canceled at 
ull points ahead of the leading edge. Cancellation of these pressure 
differences will, in turn, induce additional corrections to the perturba — 
tion velocities in region III and all portions of the wing aft of that 
region. Again, pressure differences will be introduced in the wake region 
and new corrections will have to be made which will affect the loading in 
region IV and all more rearward portions of the wing. For the present no 
tip location will be specified, except that it is, of course, outboard 
of the intersection of the first trailing-edge Mach line with the leading 
edge. The tip region will be considered after the result of the trailing- 
edge - leading-edge interference effect is determined. 

Equation (3) as originally prescribed in reference 1 (equation (51))- is 
incorrect . An errata sheet' has been issued to correct this error. 
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The initial distribution of lift to be canceled ahead of the wing (in 
region II extended) is the result of superposing an infinite-number of 
conical flow fields originating at the various points on the trailing edge. 
To cancel this distribution of lift by the superposition of conical flows 
requires that each of the lift distributions corresponding to a single 
value of a in equation (4), as well as the lift due to (Au) Q (equa- 
tion (3))> be canceled separately by an infinite number of conical fields. 

The single conical field of (Au ) 0 will be considered first. The 
velocity field to be superimposed ahead of the leading edge to cancel the 
velocities (Au ) 0 induced in the plane of the wing by the sy mme trical 
solution (equation ( 3 )) can be built up, as shown in the sketch on the left, 



of overlapping constant-velocity sectors having one edge along the leading 
edge of the wing and one along the extended ray t Q from the apex of the 
trailing edge. The magnitude of the constant velocity on each element- is 
[d(Au) 0 /dt 0 ] dt Q or, from equation ( 3 ), 



Similarly, for each oblique wake or trailing— edge element a, (see 
right-hand sketch) , a canceling field can be built up ahead of the leading 
edge by the superposition of sectors bounded by the leading edge and by 
rays t a from the apex x a ,y a of the element a and having a constant 
velocity of magnitude 
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d(Au) a 1 d (l^)(t a -m fc )-(m t -a)(l-t a ) 

— r dt„ u fi -r — cos A — dt a 

^t a n ( 1 -in t)(t a -a) 

(from equation (4)). 

Elementary Solution for the Eegion Ahead of the Leading Edge 

In general, the elementary solution required for the cancellation 
of pressure in the plane of the wing ahead of the leading edge is one 
that: 


1. Provides constant streamwise velocity over an infinite .sector 
bounded on one side by the leading edge of the w ing (extended) nrvi on 
the other by an arbitrary ray extending outward into the stream from 
some point x^y-^ on the leading edge. (See fig. 2.) 

2. Induces no vertical velocity, or downwash, on the wing. 

3. Induces no lift except on the wing and within the sector 
described in condition (l). 

At first glance these conditions would appear to be satisfied by the 
oblique solutions used at the trailing edge in the previous work, if 
properly oriented with respect to the wing, and one might expect the same 1 
form of solution to apply. In reference 1, however, it was pointed out 
that the downwash connected with this solution remained constant over the 
wing, only provided that the wing area did not include the line 
y = constant extending downstream from the apex of the element. In the 
case of the leading-edge element this condition is violated (fig.' 2) nnri 
an additional term is needed to bring the downwash to zero throughout the 
area of the wing affected by the element. 

The solution applicable to this case has been given by Lagerstrom 
(reference 2). The u component of the velocity in the plane of the w ing 
is as. follows: 



cog-1 ( ( 1+t b H m ~ t b l [ 1+t a ) 
(l+m)(t b -t a ) 


2m 

( l+m)t a 


t a ~ m) ( l+t a ) 


/ 


1+tb 
m-tb . 


( 6 ) 


In equation (6), Ub .is the constant streamwise perturbation 
velocity over the element, and tb is p times the inclination, with 
respect to the stream, of the ray from the apex x b ,y b of the element 
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through the point x,y at which the pressure is being computed. Thus, 


*b = 


3(y w y-b) 

x-x b 


(7) 


When the correction is being made for the symmetrical t railing-edge ele- 
ment only, t a is replaced by t D and Ub = [d(Au) 0 /<it 0 ] dt Q . 

For brevity, the two parts of the correction function will be 
referred to as 


c(t a ) 


r.p. 


„__i ( t a“ ffl )( 1+t bH m - t b)( 1+t a) 
cos . ■ 

(l+m)(tb-t a ) 


( 8 ) 


and 


R(t a> • ^ (iS; "/ (t a -m)(l+t a ) yg (9) 

These functions and the induced velocity (equation (6)) are plotted 
against t*^ in figure 2 for typical values of the other parameters.. 


Leading-Edge Corrections to the Local Lift 

Applying equation (6) to the cancellation of the symmetrical-wake— 
correction velocities (Au) 0 ahead of the wing results in the following 
induced— velocity increment at any point x,y on the wing: 

(Au b ) c = i r [c(t 0 ) + E(t 0 )3 dt 3 (XO) 

T o 


where t q is that value of t Q for which t-^ = — 1, and designates the 
most rearward leading-edge element containing the point x,y within its 
Mach cone. In terms of x and y, 


• = m(x+3y) 

° (x+|3y)— (l+m)c 0 


( 11 ) 


d(Au) 0 

Integration of — ^ C(t Q )dt 0 is not possible by elementary means. 

For values of m equal to or greater than O’. 4, a series expansion can be 
used for (Au) 0 (Appendix B) and a satisfactory expression found for the 
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first term of the integral (equation (10)). Calculations for m = 0.2 
indicate that integration by series may not be sufficiently accurate 
when m is small. In such a case, the first part of the integration 
may be performed graphically. ' 


The second part of equation (10) can be integrated in closed form as 
follows : 


if 


1 d ( Au )o , . 

B(t 0 ) at 0 = 


— 4m 3//£ u, 


dt f 


/ 

*(hh- ) v 


x+Py 


jtmt(l+m)K , (mt) v mx-Py |_ 


KE(k,i|r)-EF(k,\|r) 


( 12 ) 


where K and E are the complete elliptic integrals of the first and 
second kind, respectively, of the modulus 


= / 2mt(l— T n ) 

* (l-m t )(T 0 +m t ) 


(13) 


and F( k,\|r ) and E(k,ty) are the corresponding incomplete integrals, with 
the argument 


\|r = sin x 


/ T o +m t 

V 


(l4) 


The corresponding procedure for each oblique trailing-edge element 
a results in an expression for [d(Au)t,/da] da similar to equation (10). 
Again, the first term cannot be integrated in closed form; the second term 
becomes 


1 

jt 


/^•d(Au) a 


»(t a ) *t a 


where 


-4m 1— e, / ( m^-a ) ( x+Py )-( 1+m) c Q 

jt 2 (l+m)a v mx-0y 



( l + a ) ( T fl — a ) 
( l-a)(nu-a) 






(15) 


m( n^-a ) ( x+Py )-em t c D ( 1+m) 
( m^-a) ( x+py )-m t c 0 ( 1+m) 
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k = 


r (l +mt )(l-T a ) 

(l-^JCl+Ta) 


= sin" 


'-‘Til 




0(l+mt) 


(17) 

(18) 


t n = sin-i 


!(l +m t) 


(19) 


In region III (fig. l), the total leading-edge correction to the per- 
turbation velocity at any point x,y is 

(Au) b = (Au b ) 0 + f ' da (20) 

where a 0 * corresponds to the most rearward trailing-edge element the 
velocity field of which, reflected in the leading edge, can affect the 
point x,y. 


In region IV, the cancellation of the velocity increment (Au) b at 
the trailing edge must he taken into account. To do this in any rigorous 
way is obviously impractical. However, the velocities to be canceled are 
small, and it is possible to estimate the effect of their cancellation in 
many cases. Beyond region IV the velocity distribution cannot be deter- 
mined in practice, except very .approximately, by the conical— flows method. 


Calculated Results 

Load distributions have been calculated by the foregoing method for 
three combinations of taper, sweep, and Mach number as follows: 


Untapered Tapered 

m = 0.2 0.4 0.4 

nL(. = 0.2 0.4 0.6 

These values of m and m^. represent, by virtue of the Prandtl-Glauert 
transformation, a variety of practical sweep angles at Mach numbers 
between 1 and 2; as for example, 0.2 would be the value of m for a wing 
with 63° sweep of the leading edge at a Mach number of 1.075, or 75° sweep 
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at a Mach number of I.25. Similarly, m = 0.4 would correspond to 45° 
of sweep at M = 1.08, 60° at M = 1.22, or 75° at M = 1.80. The • 
trailing-edge sweep angles at these latter Mach numbers, if mt = 0.6, 
are 34°13* , 49°, and 68°, respectively. 

In figure 3> the m = 0.4 wings are shown at M =a/"2, so that 
A = cot -1 0.4. Keeping the same angle of sweep, the m = 0.2 case is 
depicted (fig. 4) as representing the same wing at a lower Mach number 
(3 = l/2). Figures 3 and 4 show also the spanwise stations at which the 
load distributions were calculated. 

The results of the calculations are presented (figs. 5, 6 , and 7) in 
the form of values of 3 (Ap/qa) where Ap is the local' lift, q the 
free— stream dynamic pressure, and a the angle of attack. In linear 
theory this quantity is connected with the streamwise component of the 
upper-surface perturbation velocity u by the relation 


= 4 — 

qa Va 

The various components of the lift, as calculated by the superposition 
method, are shown individually. 

Figure 5 presents the resulting lift distributions at the two 
stations of the tapered wing shown in figure 3* Section A— A contains the 
intersection of the trailing-edge Mach line with the leading edge, so that 
the leading-edge correction at the leading edge is zero. At points 
farther back along the leading edge, as at 3y/co =0.8, the correction 
is minus infinity. However, it is seen to increase to a small positive 
value within a fraction of the chord length at 3y/c 0 =0.8. 

At both stations it is necessary to estimate the effect of cancella- 
tion of the leading-edge correction at the trailing edge to satisfy the 
Kutta condition. Cancellation would be effected by means of oblique ele- 
ments of the type used previously (equation (4)) in canceling lift at the 
trailing edge, so that the induced lift on the wing has the same general 
shape as the primary oblique trailing-edge correction. It therefore may 
be presumed to fall along a modified inverse-cosine curve to zero at the 
boundary of the region affected. Locating this boundary, the twice- 
reflected Mach line from the trailing-edge apex, we are able to draw in a 
satisfactory estimate (dotted curve) of the correction to bring the 
pressure once more to zero at the trailing edge. 

The untapered wing with the same sweep (m = 0.4) relative to the 
Mach lines is shown in figure 6 with the load distributions calculated at 
the same stations. 

Four section lift distributions are presented (fig. 7) for m = 0.2. 
At 3y/c c = 0.15 only the rear 60 percent is influenced by the subsonic 
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trailing edge. The reflection of this influence in the leading edge 
alters the pressure over the rear 40 percent of the section. At section 
B-B the leading- and trailing-edge interaction affects the entire 
section. A further reflection of this effect in the trailing edge must 
he estimated. 

At section C— C the reflection of the leading-edge correction in the 
trailing edge covers the whole extent of the chord and any estimate of 
its magnitude would he necessarily arbitrary. Also, a second pair of ' 
reflections must be taken into account. The final pressure distribution 
has therefore been drawn as a band within which the true curve may be 
shown to lie. Its height is the error introduced at the trailing edge 
by the first leading-edge correction, except very near the leading edge, 
where an infinite negative correction is known to be introduced by the 
second leading-edge correction. The calculations were also carried out 
for 3y/c 0 = 0.45. The .margin of uncertainty was found not to have 
increased by. any appreciable amount. (See fig. 7(d).) 

It may be seen from the figures that, except in the vicinity of the 
leading edge, the magnitude of the leading-edge correction is small. The 
part due to reflecting the oblique trailing-edge correction is roughly 
the same fraction of the total leading— edge correction as the oblique 
trailing-edge correction is, at the leading edge, of the total trailing- 
edge correction. This fraction is less than 10 percent in the examples 
calculated and the resulting contribution to the leading-edge correction 
is negligible. The labor involved in computing the loading can in such cases 
be considerably reduced by omitting the calculation of the effect of 
reflecting the oblique trailing-edge correction in the leading edge. 


APPLICATION OF TWO-DIMENSIONAL FORMULAS 
TO CALCULATION OF LOAD DISTRIBUTION 

Load on Midspan Portion of Wing 

It is apparent from the calculated results that, whenever the plan 
form and the Mach number are such that the trailing-edge Mach line inter- 
sects the leading edge, the load distribution behind the Mach lines from 
the point of intersection (regions III and beyond, fig. l) resembles in 
shape the well-known incompressible load distribution over an infinitely 
long flat plate. However, the quantitative agreement does not appear 
good, particularly in the case of the tapered wing, until the load distri- 
butions over sect ions normal to the stream are examined, when a simple 
proportionality of the curves is observed. In order to determine the 
factor of proportionality it is noted that, if the curves are similar, 
each can be characterized by the strength of the s ingula rity at the 
leading edge. An explicit expression for this quantity can be obtained 
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for regions III and IV 2 of the swept— back wing, and the ratio of this 
quantity to the corresponding quantity in the two-dimensional flow is the 
required correction factor to bring the two curves into agreement. 

The leading-edge singularity in the loading on the swept— back wing 
is initially that in the triangular-wing loading. Introduction of the 
leading-edge corrections to the load in region III reduces the strength 
of the singularity there through the terms B(t 0 ) and E(t a ). (The inverse- 
cosine function is always finite.) 

It is convenient to define the strength of the leading-edge singular- 
ity as the coefficient of (mx—Py) 1 / 2 in the velocity distribution. In 
the region ahead of the trailing-edge Mach line, this coefficient is, 
from equation (l). 



reducing to 



at the leading edge. 

From equations (12) and (15)> decrements to this coefficient may be 
derived for the portion of the leading edge behind the intersection with 
the trailing-edge Mach line, as follows: 


(ACL = 


— 4muo 
itm^K^m^.) 


KE(k,t) ~ EF(k,t) 


and, for each value of a from 0 to that value a 0 which makes r a 
equal to one. 



EFO^) -/i ["keO^o) -EFO^) 


‘This expression will be shown later to hold approximately for some 


distance outboard of its theoretical limit of applicability, that is, 
beyond region IV (fig. l). * 
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where t q and 't. reduce to 


and 


T 


O 


mx 

x-c Q 


( rn^-a ) mx— n^. c Q a 
(m t -«L)x-Bi t c 0 


and the arguments and moduli of the elliptic integrals follow as for 
equations (12) and ( 15 )* 

— „ 1 J ^ 

The coefficient of (mx-(3y) in region III is therefore 


with 


Ca + (AC)o 


+ 



c Q -(i-m)x 

a ° n|.c 0 -(l-m)x 


(24) 


For convenience, we may define the nondimens ional coefficient 



In figure 8, the quantity Mx) is plotted against x/c Q for the 
wings previously discussed. Additional terms enter into the coefficient 
at the values of x/co indicated by cross-markings on the curves, where 
successive reflections of the Mach lines off the trailing edge intersect 
the leading edge. Evaluation of these terms has not been attempted; 
however, the asymptotes of the curves for the untapered wings are known 
from simple sweep theory, and it is apparent from the figure that the 
magnitude of the omitted decrements, within the practical range of aspect 
ratios, must be small. The curve for the tapered wing, without the addi- 
tional corrections, goes to zero where the wing tapers to a point, at 
x = 3c 0 . However, the complete curve, with the effect of the successive 
reflections included, would have a finite value at that point. It can 
be shown that the successive corrections enter with zero derivatives of 
successively higher order so that the error caused by neglecting oners' 
initially small. The curve in figure 8 is believed valid to the extent 
of the solid line. 
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The subsonic perturbation velocity on an infinitely long flat plate 
moving normal to its surface with circulation r is 


_T_ 

sc 



(26) 


where tj is the distance to the 
leading edge, expressed as a 
fraction of the chord c. If 
the section of the plate coin- 
cides with a section of the swept— 
back wing taken perpendicular to 
the stream (x constant), as in the 
sketch, then the chord length is 

| [mx-mt (x-c 0 )] (27) 

and 

mx-0y 

^ mx-m^x— c Q ) 

Substituting for T) in equa- 
tion (2 6), we find that the 
coefficient of (mx-£y) — 1 / e is 




i J Py-^U-Co) (29) 

which reduces at the leading edge to 

mx-mtU-Co) 

Equating this coefficient to the coefficient for the swept— back wing gives 

a_ 

T_ = 1 

nc J mx-Hn t (x-c 0 ) 

Then the two-dimensional approximation to the load coefficient on the 
outboard portions of a swept— back, wing is given by 


+ (AC) 0 + 


r dAC 


da 


( 30 ) 
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or 


u 

Va 


cr(x) 



[3y^m fc (^c 0 )]c ft 

[ mx-m^x— c Q ) ] ( mx-Py) 


( 31 ) 


— » W(x) (32 

<ia. [ mx— fflt( x— c Q ) ] ( mx-Py ) 

The closeness with which the foregoing procedure predicts the theo- 
retical loading over swept-back wings is indicated by figures 9, 10, and 
11, where the previously calculated load distributions are compared with 
those calculated by equation (32). Even in the case of the highly 
tapered wing, the agreement is seen to be good. 


At the most inboard section of the m = 0.2 wing (fig. 7(a)) there 
is, of course, no agreement over the forward portion where the flow is 
essentially conical, but behind the reflection of the trailing-edge Mach 
line in the leading edge (at 60— percent chord) the adjusted two- 
dimensional theory appears to hold fairly well. At the next section, 

B-B, the agreement is very good. At sections C-C and D-D, the curve 
calculated by the two-dimensional theory lies within the band prescribed 
for the loading by the conical— flow calculations. Since the discrepancy 
between the two-dimensional loading and the correct theoretical distribu- 
tion is already, at section B-B, less than the width of the bands in 
figures 7(c) and (d), and must diminish to zero at infinity, it may be 
supposed that the corrected two-dimensional curve is at least as satis- 
factory an approximation to the correct curves at sections C-C and 
D-D as at section B-B, and more satisfactory than can be obtained by 
the conical— flows method. 


The load distributions derived by simple sweep theory are included 
(figs. 9(l>)j 10(h), and 11(d)) to show the magnitude of the plan— form 
effect and also, in the case of the untapered wings, the curves which the 
load distributions are approaching with increasing distance from the 
plane of symmetry. 

In figures 10(b) and 11(b), comparison is also made with results 
of the s lender-wing theory of reference 4 . In the second case 
(fig. 11(b)), the Mach line from the trailing edge of the root chord, in 
the modified plan form of reference 4 , crosses the section at 12.6 
percent of the chord, causing the cusp which appears in the load distri- 
bution calculated for that wing. It is interesting to note that the 
loadings calculated by the two methods agree fairly well even when the 
wing violates the slender plan-form assumption (m=0.4). 
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Calculation of Tip Effect 

In the conical-flows method, the effect on the load distribution of 
the presence of the wing tip is calculated by canceling, outboard of the 
tip, each of the various components of the lifting pressure previously 
calculated. To do this exactly when the trailing-edge Mach line and 

leading edge intersect involves a prohibitive am ount of work. Approxi 

mate procedures have been developed (Appendix C), but the calculations 
remain tedious, while the accuracy is in some respects unsatisfactory. 

A simpler approximation, based on the assumption of cylindrical flow 
beyond the tip location, was therefore attempted. The formulas are 
developed for the tips alined with the stream, but the fundamental 
procedure may be applied to raked tips as well. 

It is assumed that the velocity distribution to be canceled in the 
plane of the wing outboard of the tip is an extension of the velocity 
distribution calculated for the tip section along lines parallel to the 
leading edge. For this purpose the approximate load distribution given 
by equation (32) is used, still further simplified by assuming a to 
remain constant at its value at the leading, edge of the tip section. The 
error involved in this assumption may be estimated by referring to 
figure 8. The maximum range of abscissa to be used is 1.0; the change in 
ordinate within that range - where the influence of the wake is effective - 
is less than 10 percent. (Where the wing is tapering to a point and a 
is changing more rapidly, the tip region is so small that the entire cal- 
culation of tip effects would probably be omitted.) 

The assumption of constant a results in a failure to cancel exactly 
the lift along the tip. The assumption of cylindrical flow, while reason^ 
able for the untapered wing (compare figs . 10(a) and (b), for example) 
would appear to be too drastic for the tapered wing, where neither the 
chord nor the loading remains constant. However, it has been shown in 
reference 1 that the major part of the tip effect results from the cancel- 
lation of the infinite pressure along the leading edge, and this part will 
be accurately calculated. The effect of the residual lift on the rearward 
portion of the tip section and in the stream will be shown, by numerical 
examples, to be small. 

The distribution of perturbation velocity at the tip station y = s 
with the simplification of constant a is, from equation (3l)> approxi- 
mately 


u(x c ,s) 


= o r 



[ps-m t (x c -c 0 )] c Q 
[mx c -mt(x c --xo) ] (mx c -ps) 


( 33 ) 
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where x c ,s are the coordinates 
of a point on the tip and a s is 
the value of cr at x = 0s/m. 

This expression may he more 
conveniently written in terms of 
the variable 




Xc-Cfts/m) 

c t 


( 34 ) 


which is the ratio of the distance 
of x c ,s back of the leading edge 
(see sketch) to the tip chord ct» 
Since 


°f ( c 'o + ^)-.V <35) 

equation (33) may be written 




If the velocity distribution x 

u is assumed to be constant beyond 
y = s along lines parallel to the 

leading edge, it may be canceled by the superposition of conical constant- 
velocity elements, having one edge along, the tip and, the other parallel to 
the leading edge, with apexes displaced along the tip by increments in 
5 C . The velocity induced at a point x,y by each such element would be 
(reference 2) 


c 

— cos 

JT 


— 1 


m+t c +2mt c 
t c — m 


( 37 ) 


where 


t c = 


P(ys) 


x— Xf 


( 38 ) 


and u c is the velocity on each sector. 

Following the procedure used in deriving equation (36) in reference 1, 
the corresponding equation may be written for the velocities induced by 
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canceling the assumed cylindrical loading: 


Au — Vm( 1 +m) ( x— x Q ) r> x ° 

Va it ' Jgs 

m 


u(x c ,s) dx r 

Va . [(l+m)x- X q-dEq ]V x q -x c 


(39) 


where x Q , s is the intersection of the Mach forecone from x,y with the 
tip. 


If the distances of x,y and x Q ,s 
ured as fractions of the tip chord, are 


hack of the leading edge, meas— 


and 



it can he shown that 


m 




(l+m)(x-x c ) = mc t (£-! 0 ) 


m 


from which equation (39) can he written (with the substitution for 
u(x c ,s) from equation (36)) 



where 


X -£t, 

, c o 


the taper ratio. 


In integrating equation (^3), three cases must he disti; 
(l) £ <1 (always true for the untapered wing), (2) 1 < | 


ished: 

^ (when the 

point x,y lies behind the trailing edge of the fictitious untapered wipg 
through the tip chord), and (3) i > (a possible condition for some 


— m 


% 

points near the trailing edge of a highly swept or tapered wing). 
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In the first case (|< 1 ) 




[KE( k‘ - (E-E^k*,^)] 


vhere K and E are the complete elliptic integrals of modulus 


1. / "• 


and E(k’,i|r 1 ) and F(k*,ijfi) are the incomplete integrals, with the 
plementary modulus 


k . . 

A-2C2 e 

v mj. »o 


and argument 



(44a) 


com- 
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In the second case ( 1 < 5 < 


C 1 


■JUfc-N 


Au 


Va 


x,y) 


~^s_ 

rtV mX. 



K - 



[KE(k,^ 2 )-EF(k,t2)] 


(44b) 


where 



In the third case 



Au, 


Va 


^y) 



[KE(kS+ 3 HM)F(k»,> a )] 


(44c) 


where 



Along the Mach line from the leading-edge tip, all three equations 
reduce to the value 


Au* _ — a s 
V a ,/ mX| 


( 45 ) 
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By the procedure just described, approximate cancellation of all 
pressure differences outboard of the tip has been effected, but the 
pressures i nd uced by such cancellation now violate the condition of zero 
lift in the wake. Cancellation of the induced pressure differences in 
the wake can be accomplished in an approximate manner, as in reference 1, 
by making use of the known value of the induced pressure at the trailing 
edge of the wing, but assuming the entire error to originate at the tip 
leading edge . 


On the trailing edge of an 
untapered wing, | = 1 and 




( 46 ) 


There is no corresponding simplifi- 
cation for the tapered wing. 

The constant— velocity elements 
to be superposed in cancellation of 
this tip-induced velocity will be 
defined by rays of slope t m/3 
from the leading-edge tip. (See 
sketch. ) Their apexes will be at 
the intersection xt,yt of the rays 
with the trailing edge, and they 
will be bounded by the rays t m on 
one side and the trailing edge on 
the other. The initial sector, 
bounded by the tip Mach line t m = — 1 
will have the constant velocity Au* 
given by equation ( 45 ). The total 
induced correction to the velocity at 
any point x,y will be 


Aeu(x,y) 


1 . 
K 




(l-m fc )(t*+l) 



1 


_ n ( 1 - t m )( t - m tHmr t m )(H j( 

5 (1-mtKt-tJ dt " 



where, if x* ,y* are the coordinates of the intersection of the tip Mach 
line with the trailing edge. 



NACA TN 1991 


21 


t* = (w) 

x— x* 

and 

t -5<22i> (49) 

x x t 

The integral in equation (4-7) is most conveniently evaluated by plotting 
the inverse cosine against Au, obtained from equation (46) for the 
untapered wing, or equations (44b) or (44c) for the tapered wing, for the 
points x-t,yt on the trailing edge. The two functions are related 
through t m , which is defined by 


P(yt-s) 
' m ” OF 

“ m 


(50) 


It is difficult to assess the accuracy of the formulas proposed for 
calculating the tip losses in the loading since no correct values are 
available for comparison. In figure 12 two streamwise section load dis- 
tributions calculated by equations (44) and (47) are compared with those 
obtained by the simplified conical— flows methods described in Appendix C. 

The latter calculations are correct along the tip Mach lines, but are 
known to fall below the correct theoretical value at more rearward points. 
The tip-induced velocity to be canceled at the trailing edge c ann ot be 
calculated with significant accuracy by the approximate method employed, 
so that no estimate could be made of the load curve within the region 
affected by such cancellation. The cylindrical— flow approximation probably 
represents as good an estimate as can be obtained at present of the theo- 
retical load distribution in that region. 


CONCLUDING REMARKS 

Beyond the spanwise station at which the Mach lines from the trailing- 
edge apex intersect the leading edge of a swept-back wing, the section 
loading takes on the general form of the two-dimensional subsonic load dis- 
tribution. A good approximation to the loading in this region can be 
obtained by applying an easily calculated correction factor to the two- 
dimensional flat— plate distribution. 

The two-dimensional approximation appears to be in good agreement 
with the theoretically more accurate loading calculated by the conical- 
flows method for tapered, as well as untapered, wings. In the case of 
tapered wings, the proposed method is able to predict the increased 
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loading in the outboard regions (similar to the effect found in subsonic 
flow) . 


The same concept of utilizing two-dimensional loading can be employed 
to obtain an estimate of the loss of loading in the zone of influence of 
the wing tip, when the tip lies wholly within the trail ing— e dge Mach cone. 
This simplification is especially valuable because even approximate appli- 
cation of the conical— flows method is tedious, and physical conditions of 
flow in the tip region are so different from the conditions assumed in the 
theory that a rough indication of the nature of the theoretical load 
should suffice at the present time as a basis for further investigations. 

The approximations of the conical— flows method, as it has been 
applied herein, should be noted: (l) In canceling lift behind either wing 
half, a small amount of upwash introduced on the other half of the wing — 
in violation of the boundary conditions — is ignored and (2), if the 
t ran ling-edge Mach line intersects the leading edge, its reflection will 
intersect the trailin g edge and, within the Mach cone from the point of 
intersection, only an estimate of the loading is feasible. If the next 
reflection (in the trailing edge) results in a second intersection of the 
Mach line with the leading edge, the leading-edge correction factor for 
adjusting the two-dimensional load distribution will also be in error. 
Calculations made to check the first of these approximations showed the 
induced-downwash angle to be less than half of 1 percent of the angle of 
attack on the wings considered. The magnitude of the inaccuracy intro- 
duced by the second limitation has been indicated by the examples pre- 
sented herein; the result is to limit the aspect— ratio Mach number com- 
binations (as expressed by the reduced aspect ratio pA) for which the 
conical— flows method will give accurate theoretical solutions. The 
corresponding error in the leading-edge correction coefficient appears, 
however, to be insignificant and it is felt that, with the exception of 
wings tapered to a point, the upper range of reduced aspect ratio PA is 
adequately covered by the 'proposed method of correcting the two- 
dimensional load distribution. There remains only the case in which the 
reduced span (p times the geometric span) is less than the tip chord. 

In this case, the interference of the tip flows with each other will, in 
conjunction with a subsonic trailing edge, create a problem not readily 
solvable by either the conical— flows method or the two-dimensional approx- 
imation. This problem falls within the scope of the slender-ving theory 
of reference 4. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Sept. 16, 1949. 
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APPENDIX A 
SYMBOLS 
General 

free— stream velocity 
free— stream Mach number 

V M 2 — 1 

density of air 
dynamic pressure ^ ^pV 2 ^ 

pressure difference between upper and lower surfaces, or local 
lift 

angle of attack, radians 

circulation on infinite flat plate 

chord of infinite flat plate (See sketch, p. 130 

Wing Dimensions 

root chord 
tip chord 
semispan 

angle of sweep of the leading edge 
taper ratio (c t /c Q ) 

Coordinates 

Cartesian coordinates in the stream direction and across the 
stream in the plane of the wing 

coordinates of apex of oblique tra iling— edge element 

coordinates of apex of leading-edge element 

coordinates of point on trailing edge, within tip Mach cone 
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x c streamy! se coordinate of apex of tip element 

x Q ■ largest value of x c at which cancellation of pressure can 
affect point at which pressure is being calculated 

I streamwise distance of x,y back from leading edge, as a fraction 

of the tip chord (equation (40)) 

! c distance of Xq behind leading-edge tip, as a fraction of the 

tip chord (equation 34) ) 

£ 0 distance of x 0 behind leading-edge tip, as a fraction of the 

tip chord (equation (4l)) 

T) spanwise distance of x,y from leading edge, as a fraction of 

the chord through x,y measured perpendicular to the stream 
(equation (28) ) 


In the following, all slopes axe measured counterclockwise from a 
line extending downstream from the apex of the wing or of the pertinent 
elementary sector: 


m 


“t 




V 

tc 


elope of leading edge = a cot A 
slope of Mach lines 

slope of trailing edge 
slope of Mach lines 

slope of ray from the origin ft y 
slope of Mach lines x 

slope of ray from trailing— edge apex 
slope of Mach lines 


3 


7 

x-c„ 


slope of ray from leading-edge tip _ ^ y- s 

slope of Mach lines x— (3s/m) 

slope of ray from apex of element a _ y-y a 
slope of Mach lines x — x 

slope of ray from apex of leading-edge element _ 

slope of Mach lines X — x^ 

% for point at which pressure is being calculated 

slope of ray from apex of tip elemen t _ Q y-s 
slope of Mach lines x— x 


n 



NACA TN 1991 


25 


a Q the value of a corresponding to a trai ling-edge element of which 

the apex lies on the Mach forecone of the point at which the 
load is being calculated 

I 

a Q * the largest value of a along which cancellation of pressure in 
the wake can affect the point at which the load is being cal- 
culated through reflection in the leading edge 

t q the smallest value of t Q along which cancellation of pressure 

ahead of the leading edge can affect the point at which pressure 
is being calculated (equation (ll)) 

T a the smallest values of t a along which cancellation of pressure 

ahead of the leading edge can affect the point at which pressure 
is being calculated (equation (l6)) 

Streamwise Components of Perturbation Velocity 

u^ .basic (uncorrected) perturbation velocity as given by solution 
for triangular wing (equation (l)) 

Uq value of u^ at a=0 (equation (2)) 

Ug constant perturbation velocity on canceling (oblique) sector 

in wake 

Ub constant perturbation velocity on canceling sector ahead of wing 

u c constant perturbation velocity on canceling sector outboard of tip 

(Au) 0 symmetrical trailing-edge correction to u A (equation (3)) 

(Au) a correction to u^ due to single oblique trailing-edge element 

(au)^ leading— edge correction to 

(AUiJq leading— edge correction to u A due to reflection of (Au) n 
(equation (10)) 

(Au) c correction to u^ due to reflection of leading-edge correction 
in wing tip (terms in E only) • 

Z^u perturbation velocity induced by canceling tip effect at the 
trailing edge (equation (47)) 
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Arbitrary Mathematical Symbols 


r.p. real part 

R radical term of leading-edge correction function (equation (9)) 

C inverse— cosine term of leading-edge correction function 

(equation (8)) 

C A value of coefficient of (mx— (3y) 1//2 in u^ at the leading 

edge (strength of the leading— edge singularity in u) 

(AC) 0 decrement in the strength of the leading-edge singularity in 
u due to reflection of (Au) 0 at leading edge 

da decrement in the strength of the leading-edge singularity in 
a u due to reflection of (Au) a at leading edge 

a leading-edge correction coefficient defined by equation (25) 

cr s value of a at leading-edge tip 


Elliptic Integrals (See references 5 and 6.) 

k modulus of elliptic integral (defined where used) 

k* the complementary modulus ( */ 1 — k 2 ) 

cp or \|r argument of elliptic integrals (defined where used) 

F(k,cp) incomplete elliptic integral of the first kind of modulus k 
and argument cp 

K,K(k) complete elliptic integral of the first kind of modulus k 
(K = F(k,§)) 

E(k,cp) incomplete elliptic integral of the second kind 

E,E(k) complete elliptic integral of the second kind of modulus k 
[E=E(k,|)J * 

K*(k) K of the complementary modulus [K'(k) = K( k’ ) ] 
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APPENDIX B 

EVALUATION OF T* — U - ° C(t 0 )dt 0 BY SERIES EXPANSION 
J T 0 dt o 

For integration of the first term of equation (10), the product 
[d(Au) 0 /dto] C(t Q ) is integrated by parts to give 



1 d(Au) 

dt 0 

o 


o 


C(t 0 )dt 0 - (Au)q 



(Au ) 0 


ac(t Q ) 
' dtg 


dto 


(Bl) 


The first term goes to zero at either limit. In the second term, tb 
must be expressed in terms of x, y, and t Q before taking the deriva- 
tive of C(t 0 ). If the constant (for any selected point) 


*o* 


Py 

X— C 0 


is defined, this derivative may be written 


dC(t 0 ) _ (l+t 0 *) >J (Tp-m) (mx— )3y) 

dt o (tp— to*) '/mc Q ( 1 +m) ( l+t Q ) ( to-^To) 


(B2) 


If this expression and the expression for (Au ) 0 given by equation ( 3 ) 
are substituted in equation (Bl), the result 


f 


d(au)o C(t„)dto - Uo(W) 


dtf 


/ 


( To-m) (mx-py) 


K( «/ l-m t 2 ) v mco(l+m) 


r 

J T 0 


F ( Vl-mt / 9) 


( tp-to* ) V ( 1 +to) ( tor“ T o) 


X 


dt r 


(B3) 


is obtained. For to close to 1 (small values of 9 ), F( Vl-mt 2 , cp) may 

be expanded into a MacLauren series in sin 9 = ( 1— 1 0 2 ) / ( 1— m^. 2 ) as 

follows : 
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From the first two terms of the series, substituted in equation • 
(B3), ve may obtain 



but the results are numerically insignificant, in most cases. 
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APPENDIX C 

. CALCULATION OF TIP LOSS IN LOADING BY CONICAL-FLOWS METHOD 

If the trai ling-edge Mach line intersects the leading edge on th 9 
wing, four components of the perturbation velocity, as calculated by 
the conical— flows method, must be canceled outboard of the tip : 3 (l) the 

basic velocity u A , (2) the trailing-edge correction (Au) Q + f da, 

J 0 

' a * 

( 3 ) the leading-edge corrections (Au) b = (Au^,) + / da > 

O J o da 

and (4) the reflection of (Au)^ in the trailing edge (in region IV, 
fig. l). This last component has not teen calculated. If it does not 
extend over more than 50 percent of the tip section, the effect of its 
cancellation at the tip can probably be neglected. 

r The third component (Au)^ can be further broken down into the 
parts arising from the radical term E in the leading— edge correction 
function and those arising from the inverse— cosine term C. The former 
components increase to infinity at the leading edge and will be treated 
in the same way as the basic velocity u^. The latter remain finite 

and exactly cancel the trailing— edge corrections (Au) n +/ da 

J da 
o 

at the leading edge of the tip section. Although these components (the 
C component of (Au)^ , and the trailing-edge correction) do not cancel 
each other elsewhere on the tip section or in the stream, the variation 
of each is fairly small in the region affecting the pressures on the 
wing, and, since the net induced effect of canceling them is zero and has 
zero slope along the tip. Mach cone, it will be assumed that they cancel 
each other completely. This assumption will result in somewhat too low 
a loading near the trailing edge. 

Cancellation of the first component, the basic velocity u^, has 
been covered in reference 1. The result is a reduction in u given by 


Au = 


^ 23 ( 3 — y) 


rm 2 ! 2 — 3 2 y 2 


£_ke' 

K* \ 2 


+ KE(k». 


If the tip lies beyond region IV (fig. l) further small components 
enter, but these will not be considered- in the present analysis. 
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(equation (37) of reference l) , where 
s is the semispan. 


k = 


' ( m-a D ) ( 1-m) 
2m(a 0 +l) 


\|r = sin 


1 


/a n (mx+3y) 


3s 


3s(a 0 +m) , a ° x+3(y-s)^ 


K* = K(k»),and E» = E(k»). 

The part of the leading-edge correc- 
tion associated with R(t 0 ) may he 
canceled by the following procedure: 

For any one value of t Q (see sketch), 
the correction to u(x,y) will be 



d ( Au )c d t 0 = - d(A ^)a d t 0 limit fR(t 0 ,t t ) cos" 1 ^ b+t c +2t ^ t c + 
dt it 2 dt t-u->m fc - ^b 


/*<=* _1 & -1 

/ cos 

Jd-_ • St. 


-1 t-u'+t„ +2t-ut. 


dt 

tc-tb ■ 1 


(C2) 


where t^ = ^ s ~* (see sketch). - 


x -" x b 


Integrating by parts results in the elimination of the first term 
so that 


d(Au) 

dt 


-‘c“ - 1 


ic at# = i, S^Blo dt 0 ^ coa- VWb dt. 


Equation (C3) can be integrated analytically, giving 


(C3) 


d(Au). 


dt^ = 


2m 


dt o ° n ( 1+m ) t c 


/ ( t 0 ~ m) ( l+t Q ) ( l+tfr* T d ( Au ) 
v- m— t * c 


(C4) 


where 


t * = 


3(y-y b ) 


X— X-, 
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To obtain the total correction to u(x,y) for the- cancellation of 
the term in R(t Q ) of the leading-edge correction, equation (c4) may be 
integrated graphically by plotting against (Au) q the remainder of the 
expression in the right — hand side. An approximate correction for the 
reflection of the terms in R(t & ) can be included by plotting against 
the values of 


dAu 


r ° 

Au = (Au) + 

0 J da. 

o 


da 


along the leading' edge. Then 

(Au) = 2m / ( tp-m) ( l+t 0 ) ( 1+t-b* ) 

c Jt(l+m) J. t« v 


1 


V 


dAu (C5) 


where t Q , t^*, and Au are all computed for specified locations of 
Xj ) ,y-j 5 on the leading edge and the extension of the leading edge to 
the point such that t-^* = — 1. 

In practice, equation (C5) appears to result in a nearly constant 
correction to the section loading, as shown in figures 12(a) and (b). 
The work may be further simplified, therefore, by performing only the 
relatively simple calculation for a point on the Mach line. 
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Figure I. - Plan view of central portion of swept- back 
wing showing pattern of Mach lines arising at 
leading and trailing edges. 
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-Uncorrected 





-Final 

r* nr* 


Leading-edge correction 


Correction for 

Kutta condition- 


Oblique trailing -edge 
■ — w correction - 


Symmetric trailing -edge correction- 


25 50 75 ~Wt 

Percent chord 

(a) Section A- A, Ay/c a -0.667. 

Figure 5. - Load distribution over two stream wise sections 
of the tapered swept -back wing calculated by the 
conical-flows method . 


NACA TN 1991 



NACA TN 1991 


39 



the untapered swept -back wing, m • 0. 4, calculated by 
the conical-flows method. 


ko 
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-Uncorrected 


Final fi AJL.{ 
Qa 



L eading - edge correction 


Oblique trailing- edge correction 


Symmetric t railing-edge correction 


\ ■ Correction for Kutta condition 

0 ^ £>d 75 Jdo 

Percent chord 

(a) Section A -A; /3y/c g - 0.15. . 

Figuse 7. - Load distributions over four streamwise sections of the 
untapered swept-back wing , m-0.2, calculated by the conical 
flows method. 
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Figure 8. - Leading-edge velocity coefficient cr(x) for correction of two- dimensional 
load distribution. 
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Figure 9- Concluded. 
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Figure // - Continued . 
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Figure 12.- Load distribution over section near tip as cal- 


culated by conical - flows method and by assuming cylindrical 
flow near tip. 
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(b) Tapered wing, /BA* 2.61. Section at /By * 0.8 c 0 
or 94-percent semispan. 

Figure !2r Concluded. 
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